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Abstract 



Aims. We have investigated a frequency-dependent shift in the absolute position of the optically thick apparent origin of parsec-scale 
jets ("core shift" effect) to probe physical conditions in ultra-compact relativistic outflows in active galactic nuclei. 
Methods. We used multi-frequency Very Long Baseline Array (VLBA) observations of 191 sources carried out in 12 epochs in 
2006 within the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) program. The observations were 
performed at 8.1, 8.4, 12.1, and 15.4 GHz. We implemented a method of determining the core shift vector based on (i) image 
registration by two-dimensional normalized cross-correlation and (ii) model-fitting the source brightness distribution to take into 
account a non-zero core component offset from the phase center. 

Results. The 15.4-8.1, 15.4-8.4, and 15.4-12.1 GHz core shift vectors are derived for 163 sources, and have median values of 128, 
125, and 88 fjas, respectively, compared to the typical measured errors of 50, 51, 35 yuas. The effect occurs predominantly along the 
jet direction, with departures smaller than 45° from the median jet position angle in over 80% of the cases. Despite the moderate ratio 
of the observed frequencies (<2), core shifts significantly different from zero (>2<x) are detected for about 55% of the sources. These 
shifts are even better aligned with the jet direction, deviating from the latter by less than 30° in over 90% of the cases. There is an 
indication that the core shift decreases with increasing redshift. Magnetic fields in the jet at a distance of 1 parsec from the central 
black hole, calculated from the obtained core shifts, are found to be systematically stronger in quasars (median B\ as 0.9 G) than those 
in BL Lacs (median B t a 0.4 G). We also constrained the absolute distance of the core from the apex of the jet at 15 GHz as well as 
the magnetic field strength in the 15 GHz core region. 

Key words, galaxies: active - galaxies: jets - quasars: general - radio continuum: galaxies 



1. Introduction 

Bipolar relativistic outflows (jets) in active galactic nuclei 
(AGN) are formed in the immediate vicinity of the supermassive 
central black hole and become detectable at distances of >100 
gravitational radii (R e - GMbhlc 2 ) at millimeter wavelength s 
(|Junoretal.iri999l: lLobanov & Zensusll2007l: lHada et al.ll20lTl) . 
The jets take away a substantial fraction of the energy and angu- 
lar momentum stored in the accr etion flow (|Huieirat et aT1l2003l) 
and s pinning central bla ck hole ( Koide et aO 2002 ; iKomissarovl 
12005b . As discussed by IVlahakis & Konigll ( 120041) . a poloidal- 
dominated magnetic field embedded in the accretion disk or in 
the black hole ergosphere is wound-up into toroidal loops that 
may provide effective jet collimation via hoop stress and accel- 
erate the flow by magnetic pressure gradient up to a distance of 
~10 3 - 10 5 R g . 

Very Long Baseline Interferometry (VLBI) observations 
provide us with the perfect zoom-in tool to explore AGN 
jets with a milliarcsecond angular resolution corresponding to 
parsec-scale linear resolution. Typically, the parsec-scale radio 
morphology of a bright AGN mani fests a one-sided jet struc- 
ture due to Doppler boosting (e.g., iBlandford & Konigll 1 19791 : 



iKellermann et alJ 120071: iLister et alJ 120091) that enhances the 
emission of the approaching jet. The apparent base of the jet 
is commonly called the "core", and it is often the brightest and 
most compact feature in VLBI images of AGN. The VLBI core 
is thought to represent the jet region, located at the distance 
'"core t° the central engine, at which its optical depth reaches 
t v « 1 at a given frequency. At short mm-wavelengths the 
core may also be the first recollimation shock downstream of 
the r = 1 surface instead of the surface itself. This does not 
affect our analysis, which uses longer wavelengths. Thus, the 
absolute position of t he radio core is frequency-dependent an d 
varies as r core oc v~ 1/k ' (IBlandford & Koniglll979r.lKonig]|l98ll) . 
i.e., it shifts upstream at higher frequencies and downstream 
at lower frequencies (the so-called "core shift" effect). The 
first core s hift measurement from VL BI observations was per- 
formed bv lMarcaide & Shapiro! dl984l). Recent multi-freq uency 



studies of the core shift effect dO' Sullivan & Gabuzda 200S 
iFromm et all 120101: ISokolovskv et al.1 1201 it lHada et alJ (201 lh 
showed that k r « 1 in most sources and epochs. This is consis- 
tent with fhe lBlandford & Konigll (fl979h model of a synchrotron 
self-absorbed conical jet in equipartition between energy densi- 
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ties of the magnetic field and the radiating particle population. 
Nonetheless, departures in k r from unity are also possible and 
can be caused by pressure and density gradients in the jet or 
by external absorption fr om the surrounding medium dLobanovl 
[l998t iKadler etaLll200l . 

The frequency-dependent offsets of the core positions can 
be used for astrophysical studies of ultra-compact AGN jets to 
calculate the magnetic fields, synchrotron luminosities, total (ki- 
netic and magnetic field) power, maximu m brightness tem pera- 
ture and geometrical properties of the jet dLobanovlll998l) . The 
core shift effect also has immediate astrometric applications. 
A typical shift between the radio (4 cm) and optical (6000 A) 
domains for distant quasars i s estimated to be at the level of 
0.1 mas dKovalev et al.1 120081) . which is comparable with the 
expected positional accuracy of the GAIA astrometric mission 
dLindegren & Perrvmanl 19961) . Thus, the core shifts are likely to 
influence not only the positional accuracy of the radio reference 
frame but also an alignment of optical and radio astrometry cata- 
logs. Moreover, it is natural to expect that opacity properties are 
variable on a time scale from months to years due to the contin- 
uous emergence of new jet components, and esp ecially during 
strong nuclear flares. Therefore, as discussed by iKovalev et alj 
d2008l) . a special coordinated program is required to perform 
multi-frequency and multi-epoch VLBI observation of a pre- 
selected source sample to investigate the problem of core shift 
variability. 

A major difficulty in measuring the core shift is the ac- 
curate registration of the VLBI images taken at different fre- 
quencies. The problem stems from the loss of absolute po- 
sition information in the standard VLBI data reduction path, 
which involves self-calibration of the station phases. Several ap- 
proaches have been presented to overcome this difficulty and 
measure core shifts. One of them is based on relative VLBI 
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This particular technique is resource-consuming and has been 
used for a limited number of s ources only. An o ther approach 
is the self-referencing metho d dLobanovl 11998; Ko valev et alj 
120081: ISokolovskv et al.l[20TTh . in which the core shift is derived 
by referencing the core position to bright optically thin jet fea- 
tures whose positions are expected to be achromatic. Although 
this method has provided the majority of known core shift mea- 
surements, it has a certain limitation. It cannot be applied for 
faint or smooth jets that lack compact bright feature(s) well sep- 
arated from the core at different frequencies. A proper align- 
ment of the optically thin parts of the jet can also be accom- 
plished by two-dimensional cross -correlation of the im ages, ini- 
tially suggested and performed by Walker et al. ( 2000) for multi- 
frequency VL BA observations of 3C 84 . The algorithm was also 
discussed by ICroke & Gabuzda J 2008b and Fromm et al. (in 
prep.), and applied by |Q' Sullivan & Gabuzdal d2009h to obtain 
core shifts in four BL Lac objects. This approach, in conjunc- 
tion with source model fitting, presents a more widely applicable 
method for deriving core shifts (see Sect. [3] for detailed discus- 
sion), which we use in this paper. Another alte r native indirect 
method recently proposed bv lKudrvavtseva et al.l d201 lb is based 
on an analysis of time lags of flares monitored with single-dish 
observations. Although it has obvious limitations on the epoch 
at which the core shift can be measured, the method is promis- 
ing for highly compact sources, which pose problems for other 
opacity study methods due to the lack of optically thin jet struc- 
ture. It is noteworthy that all of the aforementioned techniques 
provide a comparable accuracy level. To date, only two core shift 



studies dKovalev et al.ll2008b ISokolovskv et al.ll201 ll) have been 
carried out on large samples. They have shown that the effect is 
significant for many sources. 

In this paper, we measure frequency-dependent shifts in the 
absolute core positions and study the statistical properties of the 
detected core shift vectors by using a large sample of sources 
from the MOJAVE (Monitoring Of Jets in Active galacti c nu- 
clei with VLBA Experiments) program dLister et al.| [2009). We 
also analyze systematics and discuss the uncertainties of the two- 
dimensional cross-correlation technique, investigating its prop- 
erties for different jet morphologies. We constrain the basic 
physical properties of the jets, such as the magnetic field strength 
in the core region and at the true base of the flow, the distance 
from the jet apex to the radio core, as well as the estimate of the 
central black hole mass from the derived core shifts. 

Throughout the paper, we assume the power index k r = 
1, i.e., r C0le °c v _1 (see model assumptions for this case 
above). We use the ACDM cosmological mod el with Hp = 
71 km s" 1 Mpc" 1 , Q,„ = 0.27, and C2 A = 0.73 dKomatsu et al.1 
120091) . All position angles are given in degrees from north 
through east. 



2. Observations and data processing 



The MOJAVE project dLister et al.ll2009l) is a long-term VLBA 
program aimed at investigating the structure and evolution of 
extragalactic relativistic radio jets in the northern sky. The full 
monitoring list currently consists of about 300 sources, and in- 
cludes a statistically complete, flux-density limited sample of 
135 AGN, referred to as MOJAVE- 1. In addition to the pro- 
gram's usual single-frequency setup at 15.4 GHz, 12 monthly 
separated epochs of observations during 2006 were carried out 
simultaneously also at 12.1, 8.4, and 8.1 GHz. The observations 
were made in dual circular polarization mode, with a bandwidth 
of 16 MHz at two lower bands and 32 MHz at two upper bands, 
and recorded with a bit rate of 128 Mbit s -1 . In total, 191 sources 
were observed. 

The initial calibration was performed wi th the NRAO 
Astronomical Image Processing System (AIPS) dGreisenll2003l) 
following the standard techniques. All frequency bands were 
processed separ ately throughout the data reduction. C LEANing 
dHogbomlll974), phase and amplitude self-calibration dJennisonl 
195 81: iTwiss et al.lll960l) . were performed in the Caltech Difmap 
dShepherdl 119971) package. In all cases a point-source model 
was used as an initial model for the iterative procedure. Final 
maps were produced by applying natural weighting of the vis- 
ibility function. For a more detailed discus sion of the data re - 
duction and imagin g process schemes, see iLister et ail d2009t) : 
iHovatta et all d2012l) . 

The structure of each source at each frequency band was 
model-fitted in the visibility (u, v) plane in Difmap using cir- 
cular and elliptical Gaussian components. To achieve matched 
resolution in all bands in the image plane, we appropriately cut 
the long baselines from the 15.4 and 12.1 GHz interferometric 
visibility data sets and short baselines from 8.1 and 8.4 GHz 
data sets. For each source, all maps were restored with the same 
beam size taken from lowest frequency (8.1 GHz) data using a 
pixel size of 0.03 mas, and these images were cross-correlated 
to register them as explained in Sect. [3] 

This set of observation s was also used to in vestigate the jet 
Faraday rotation measures dHovatt a et al. 20 12j) and spectral in- 
dex distributions (Hovatta et al., in prep.). 
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3. Method for measuring the core shift 

Registration of two VLBI images taken at different frequencies 
provides a shift between the image phase centers, such that they 
are co-aligned to the same position on the sky. This image shift 
is not equal to the core shift we want to measure. In general, the 
presence of the jet structure we detect results in a non-zero off- 
set of the core from the phase center. As the core components 
significantly dominate in total flux, the magnitude of the offsets 
is typically small, but at the same time not negligible. When a 
distant jet component is brighter than the core, the offset can be 
large, as in the case of the quasar 0923+392, whe re the offset 
is abo ut 2.6 mas. Extreme cases are discussed in iPetrov et alJ 
d20Tlh . Moreover, these offsets are different at different frequen- 
cies for a given source, and become statistically larger at lower 
frequencies due to spectral properties of the jets. FigureQ]shows 
the core offsets at 15.4 and 8.1 GHz, with the median values 
being 36 and 8 1 //as, respectively. 

The core position departures from the phase center thus 
have to be taken into account to derive the core shifts. Because 
the image shift IS is independent of the core component posi- 
tion, it measures the vector sum of the absolute core shift CS 
(\CS\ = Ar core VlV ,) and the difference in the coordinates (offset 
shift relative to the map center) OS: 
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Figure 1. Phase center offsets of the core components at 8.1 
GHz versus that at 15.4 GHz taken from model fits. 



IS = CS + OS, (1) 

from which the magnitude and direction of the core shift vector 
can be readily calculated. 

We used the fast n ormalized cross-correlation (NCC) algo- 
rithm bv lLewisI d 1995b to register the images across the frequen- 
cies. The algorithm allows one to apply frequency domain meth- 
ods to calculate the unnormalized cross correlation and then ef- 
ficiently normalizes it by using precomputed integrals of the im- 
ages over the search area. Thus, spatial domain computation of 
the cross correlation function is not needed and for large images 
the decrease in computing time is significant. The features that 
were matched between the images were selected from the opti- 
cally thin part of the jet and assumed to have constant spectral 
index across them. The effect of possible spectral index gradients 
across the features is discussed in Sect. [5] In all cases, the image 
shifts IS, obtained with NCC, were verified by visually inspect- 
ing the corresponding spectral index images before and after 
the alignment, because the latter are extremely s ensitive to the 
in)acc urate image alignment, as was shown by iKovalev et alJ 
2008). The spectral index images are presented and discussed 
by Hovatta et al. (in prep.). The image shifts between 15.4 GHz 
and other bands are found to be within a range of and 1.11 mas, 
with a median value of 0. 1 3 mas. The extreme value of 1 . 1 1 mas 
was detected in 3C 273 between 15.4 and 8.4 GHz, where the 
peak of brightness at those two frequencies corresponds to dif- 
ferent features in the source structure. In some rare cases (e.g., 
1928+738 at epoch 2006 Apr 28 of and 2128-123 at epoch of 
2006 Oct 6), when an optically thin, bright, compact jet compo- 
nent dominates in flux density, the image shift was measured to 
be zero, as expected for achromatic jet components. The accu- 
racy of the two-dimensional cross-correlation technique is dis- 
cussed in Sect. [5] 

The magnitude of the core offset difference vector OS be- 
tween 15.4GHz and other bands ranged between and 1.09 mas, 
with a median value of 0.05 mas. The maximum value of 
1.09 mas holds also for quasar 3C 273 due to the same reason. 



4. Core shift measurement results 

Substituting the results of model fitting and two-dimentional 
cross correlation into Eq. ([1]), we calculated the magnitude 
Ar core , VlV2 and direction CS of the 15.4-8.1, 15.4-8.4, and 15.4- 
12.1 GHz core shift vectors for 160, 158, and 147 sources 
(Table [TJ, respectively. For the other 31, 33, and 44 sources we 
could not measure the respective 15.4-8.1, 15.4-8.4, and 15.4- 
12.1 core shifts, mostly due to the weakness of their jet emission 
(especially at 15.4 GHz). This made the cross-correlation tech- 
nique inapplicable, since there was no sufficiently large optically 
thin emission structure for feature matching. We also excluded 
those sources, mostly nearby galaxies listed in Table |2j whose 
core region was com plex (e.g., 3C 84, M 87; the core shift in 
M 87 was studied by lHada et alJl20lH using phase-referencing 
VLBA observations) or for which the identification of the core 
component was unclear (e.g., 0108+388, 1509+054). The max- 
imum and median magnitude of the derived 15.4-8.1, 15.4-8.4, 
and 15.4-12.1 GHz core shift vectors in angular and linear scale 
are summarized in Table [3] As seen from Table [3] the median 
values of the 15.4-8.1 and 15.4-8.4 GHz core shifts are compa- 
rable, while the 15.4-12.1 GHz ones are statistically smaller, as 
expected. In angular scale, these values are of about 8% of the 
corresponding FWHM beam size at 8.1 GHz. 

In Fig. [2] we present plots of the derived core shifts in polar 
coordinates, where the head of each core shift vector represents 
the core position at lower frequency, while all core positions at 
higher frequency are placed at the origin. The polar axis corre- 
sponds to the median jet direction 0j et calculated from position 
angles of the jet components with respect to the core component, 
using a corresponding model fit at 15.4 GHz. Thus the polar co- 
ordinates of the head of each vector represent the magnitude of 
the core shift vector, Ar core V|V ,, and the angular deviation from 
the jet direction, CS - 6>j et . The shift effect occurs predominantly 
along the jet direction. In more than 80% of cases, the core shift 
vectors deviate less than 45° from the median jet position angle. 
Statistically, the larger core shifts have better alignment with the 
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Table 1. Derived core shift vectors. 



Source 


Epoch 


Median 


15.4-8. 


1 GHz core shift 


15.4-8 


.4 GHz core shift 


15.4-12.1 GHz core shift 






jet PA 


PA 


total 


proj 


PA 


total 


proj 


PA 


total 


proj 






(deg) 


(deg) 


(mas) 


(mas) 


(deg) 


(mas) 


(mas) 


(deg) 


(mas) 


(mas) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


AAA') fMLC 

0003—066 


2006-0 /-0/ 


-82.1 


-60.3 


0.035 


0.033 


-25.8 


0.019 


0.011 


-62.5 


0.015 


0.014 


0003+380 


2006-03-09 


117.0 


77.2 


0.134 


0.103 


79.6 


0.139 


0.110 


77.2 


0.124 


0.095 


0003+380 


2006-12-01 


116.2 


115.5 


0.063 


0.063 


121.7 


0.106 


0.106 


103.0 


0.046 


0.044 


0007+106 


2006-06-15 


-67.6 


-88.5 


0.008 


0.007 


-20.1 


0.011 


0.007 


146.2 


0.008 


-0.007 


0010+405 


2006-04-05 


-31.8 


-38.3 


0.013 


0.013 


64.7 


0.008 


-0.001 








0010+405 


2006-12-01 


-32.6 


-0.3 


0.005 


0.004 


39.1 


0.005 


0.001 


-89.9 


0.010 


0.006 


0055+300 


2006-02-12 


-50.1 


-45.2 


0.179 


0.179 


-10.2 


0.083 


0.064 


-61.5 


0.053 


0.052 


0106+013 


2006-07-07 


-125.2 


-113.7 


0.005 


0.005 


-139.2 


0.005 


0.005 


20.2 


0.002 


-0.001 



Notes. (1) IAU name (B 1950.0); (2) epoch of observations; (3) 15.4 GHz median jet position angle; (4) position angle of the 15.4-8.1 GHz core 
shift vector; (5) magnitude of the 15.4-8.1 GHz core shift vector; (6) 15.4-8.1 GHz core shift vector in projection on the median position angle; 
(7), (8), (9), and (10), (11), (12) the same as (4), (5), (6) but for 15.4-8.4 GHz and 15.4-12.1 GHz core shifts, respectively. TableQ]is published in 
its entirety in the electronic version of Astronomy & Astrophysics. A portion is shown here for guidance regarding its form and content. 




Table 2. Sources excluded from the analysis due to unclear core 
position. 
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Figure 2. Polar plots of the 15.4-8.1 (top), 15.4-8.4 (middle), 
and 15.4-12. 1 GHz (bottom) core shift vectors. The polar axis is 
pointing to the right and co-aligns with the median jet direction 
of a source. Radial distance is given in mas. Dotted lines are 
drawn at intervals of 30°. 



Source 


Alias 


Opt. class" 




0108+388 




G 


0.668 


0238-084 


NGC 1052 


G 


0.005037 


0316+413 


3C 84 


G 


0.0176 


0429+415 


3C 119 


Q 


1.022 


0710+439 




G 


0.518 


1228+126 


M 87 


G 


0.00436 


1404+286 


OQ 208 


G 


0.077 


1509+054 




G 


0.084 


1607+268 


CTD93 


G 


0.473 


1957+405 


Cygnus A 


G 


0.0561 


2021+614 


OW637 


G 


0.227 



Notes. As compiled bv lLister etall d2009h 



Table 3. Core shift statistics. 



Core shift 




anj 


*ular 




linear 




N 


max 


med 


N 


max med 






(pas) 


(pas) 




(pc) (pc) 


15.4-8.1 GHz 


160 


525 


128 


154 


3.108 0.658 


15.4-8.4 GHz 


158 


449 


125 


152 


2.927 0.643 


15.4-12.1 GHz 


147 


286 


88 


142 


2.012 0.429 



jet direction because (i) they are less influenced by random er- 
rors and (ii) the core shift takes place along the jet in most cases. 
The weighted average of # cs - 6>j e t is close to zero. Significant 
angular deviations of the core shift vectors from the median jet 
direction may take place in sources with substantial jet bending, 
either within an unresolved region near the VLBI core, or in the 
outer jet, thus affecting the median jet position angle. We have 
also analyzed distributions of the angular deviation between the 
core shift vectors and (i) inner jet direction determined as a po- 
sition angle of the innermost jet component at 15.4 GHz and 
(ii) flux density-weighted average of the position angles of all 
Gaussian jet components at 15.4 GHz. In both these cases the 
scatter was larger, indicating that the median jet position angle 
is a better estimate for the direction of the outflow for the major- 
ity of sources. 

Our median 15.4-8.1 GHz core shift of 0128 mas exceeds 
that of 0.080 mas reported by ISokolovskv et al.l d2011l) . which 
was based on a smaller sample of 20 sources, for which the core 
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shifts were derived using the self -referencing method. But if the 
15.4-8. 1 GHz core shifts fo r the 20-source sample a re calculated 
from the fitted hyperbolas dSokolovskv et al.ll2.01 ll) . which pro- 
vide more accurate core shift values, the corresponding median 
yields 0. 127 mas, which agrees well with the median of our sam- 
ple. 

The sources with the largest angular shifts are all at z < 1, as 
shown in Fig.[3](top), where we plot the 15.4-8.1 GHz core shift 
against redshift values. All measurements are confined under the 
aspect line corresponding to 3.1 pc, the maximum linear 15.4- 
8. 1 GHz shift (see TableO. In contrast, in linear projected scale 
(Fig. [3] bottom), low -redshift sources are characterized by small 
shifts due to quick falling of the scale factor, while more distant 
sources have larger shifts. We have found evidence that the an- 
gular core shift decreases with increasing redshift by binning the 
data into nine equally populated bins, though we cannot claim 
the dependence to be highly significant due to the large scat- 
ter in measured core shifts, which stems from different physical 
conditions and viewing angles in different sources at the same 
redshift. The plots for 15.4-8.4 and 15.4-12.1 GHz core shifts 
are qualitatively similar to Fig. [3] 

Additionally, we tested whether the core shift measurements 
are affected by Um ited angular resolution blending following the 
approach used by iKovalev et al.l d2008l) . If present, this effect 
would preferentially increase the magnitude of the core shift vec- 
tors when they are better aligned with the major axis of the in- 
terferometric restoring beam resulting in a U-shape dependence 
between Ar coreiVlV2 and |# C s _ #bpal, where #bpa is the position angle 
of the major axis of the beam, and -90 < #b pa < 90. We found 
no such trend in the 15.4-8.1 GHz measurements, confirming 
that the registered core shifts are not dominated by blending. 



5. Accuracy of the method 

5. 1 . Random errors 

As seen from Eq. (Q}, the uncertainty of the core shift cr VlVo is 
determined by errors of the core positions and image registra- 
tion. However, individual a posteriori estimates of the image 
registration accuracy are problematic, since (1) in our case the 
values of the applied similarity criterion (NCC) are not normally 
distributed, making an estimation of random errors difficult, and 
since (2) systematic errors can only be addressed by simulations. 
Therefore, we used a statistical approach to assess the typical 
random core shift error in our sample. If we assume that (i) the 
core shift occurs along the jet and (ii) the errors are random in 
direction and comparable to each other, then the standard devi- 
ation of the transverse projections of the core shift vectors onto 
the jet direction yields the typical error. In Fig. [4] we plot the cor- 
responding 15.4-8.1, 15.4-8.4, and 15.4-12.1 GHz distributions, 
from which we findcr 15 4 _g 1GHz = 50yuas, cr 15 . 4 _ 8 .4GHz = 51juas, 
and <r 15 .4_i2.i ghz = 35 yt/as. From these values, we determine that 
in 57%, 59%, and 58% of cases the respective core shifts are sig- 
nificantly (>2<x) different from zero, more than 90% of which in 
turn deviate less than 30° from the median jet direction. The de- 
rived error estimates are conservative, since in some cases the 
angular deviations of the core shift vector from the median jet 
direction can be real, for instance, in curved jets. 

An alternative way to estimate a typical random error in 
core shift is based on the fact that 8.1 and 8.4 GHz bands are 
closely separated, but were processed independently. Therefore, 
the 15.4-8.1 and 15.4-8.4 GHz core shifts are expected to 
be virtually the same and the non-zero difference between 
them is due purely to errors. This approach was also used in 
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Figure 3. Top: 15.4-8.1 GHz core shift in angular scale versus 
redshift. The dashed lines correspond to a fixed projected linear 
scale of 0.1, 1, and 3.1 pc. The latter curve envelopes all the 
measurements. Bottom: 15.4-8.1 GHz core shifts in projected 
linear scale versus redshift. Typical error values are discussed in 
Sect. [5] 



ISokolovskv et al.l (1201 ll) . Constructing a distribution of these 
differences and calculating its standard deviation, we found 
ci5.4-8.ighz ~ 0"i5.4-8.4GHz ~ 55 //as, which is consistent with 
the error estimates obtained using the first method. 
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Figure 5. Naturally weighted total intensity CLEAN image of 
0716+714 at 15.4 GHz (top). The alignment error (in beam 
widths) along the jet as a function of frequency ratio (middle). 
The alignment error (in beam widths) perpendicular to the jet 
(bottom). The error is defined as r Vl -r v , where r Vi is the position 
in the low-frequency image and r Vl is the corresponding position 
in the high-frequency image. The different spectral index gradi- 
ents are shown in different colors: -0.1 mas -1 (black squares), 
-0.2 mas -1 (blue crosses), and -0.3 mas~' (red circles). 
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Figure 6. Naturally weighted total intensity CLEAN image of 
3C273 at 15.4 GHz (top). The alignment error (in beam widths) 
along the jet as a function of frequency ratio (middle). The align- 
ment error (in beam widths) perpendicular to the jet (bottom). 
The error is defined as r Vl — r v , where v Vl is the position in 
the low-frequency image and r Vl is the corresponding position 
in the high-frequency image. The different spectral index gradi- 
ents are shown in different colors: -0.1 mas -1 (black squares), 
-0.2 mas -1 (blue crosses), and -0.3 mas -1 (red circles). 



Since there is some freedom in selecting the jet feature to 
be matched in NCC, there is a possibility that user decisions af- 
fect the image registration results. We tested the robustness of 
the registration algorithm to "user bias" by having two different 
people separately perform the image alignment for one of the 
observing epochs. They both had similar instructions regarding 
the selection of the matched feature, i.e., (1) it should be opti- 
cally thin, (2) and it should have as much structural variation as 
possible. The distribution of the differences in IS had a standard 
deviation of w 50 //as for N - 24- pairs of images. This also 
closely matches the typical random error estimated above. 



5.2. Systematic errors 

So far we have discussed random errors, but the image registra- 
tion method is also prone to systematic errors that may bias our 
core shift measurements. Namely, the assumption of a constant 
spectral index across the matched features does not necessarily 
hold in real jets. Indeed, it is known that spectral index gradi- 
ents occur along the jet. The spectral index images between 8.1 
and 15.4 GHz typically show gradients in the optically thin part 
of the jet from -0.3 mas -1 to +0.1 mas -1 with an average of 
~ -0.1 mas -1 (Hovatta et al., in prep.). 
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We tested the systematic effect that such gradients may have 
on image registration results by performing simulations. For the 
simulations we selected five sources with different jet morpholo- 
gies: NGC315 (straight, long and very smooth jet), 0716+714 
(straight, short and smooth jet; Fig. [5] top), 3C 273 (long, wig- 
gling jet with prominent knots; Fig. |6l top), BLLac (wide, 
curved jet having one prominent knot downstream of the core), 
and 3C 454.3 (complex, curved jet with knots). For each source 
we created three sets of simulated images at six different fre- 
quencies exhibiting three different spectral index gradients along 
the jet: —0.1, -0.2, and -0.3 mas -1 . The simulated images were 
based on a real 15.4 GHz image of a given source, to which a 
constant spectral index gradient along the jet was applied and 
new images at frequencies of 1.25, 1.50, 1.75, 2.00, 2.25, and 
2.50 times the original frequency were calculated. Finally, ran- 
dom noise at the same level as in the original image was added 
to simulated images to ensure that background noise patterns 
between the images do not correlate. The original image and the 
simulated one were then registered using NCC. Note that this 
simulation setup provides a worst-case scenario in the sense that 
the gradient is assumed to be constant for the whole jet length, 
whereas in the real jets this is typically not the case. 

The simulation results show that spectral index gradients 
along the jet can indeed affect the registration results and that 
the systematic error introduced this way depends on the jet mor- 
phology, the magnitude of the gradient, and the frequency ratio 
of the registered pair of images. In 3C273 and BLLac, which 
have significant structural detail in the optically thin part of the 
jet, the errors in the registered shift along the jet are less than 5% 
and less than 3% of the beam size, respectively, for all values of 
the gradient and for frequency ratios lower than 2. For 3C273, 
a gradient of -0.1 mas -1 results in an error that is less than 1% 
of the beam size (Fig. [6j middle). The jet in 3C 454.3 also has 
significant structural detail, and the systematic errors for a gra- 
dient of -0.1 mas" 1 stay below 4% of the beam size. However, 
for steeper gradients, the errors increase significantly, being less 
than ~10% of the beam size for frequency ratios below 2. The 
featureless jets of NGC 315 and 0716+714 are the most prone 
to systematic errors caused by spectral index gradients along the 
jet: even the flattest gradient results in errors in the range of 10- 
18% of the beam size. In 0716+714 the gradients of -0.2 and 
-0.3 mas" 1 result in errors of 24% and 27% of the beam size, re- 
spectively, at a frequency ratio of 2 (Fig. [5] middle). In NGC 315 
the steeper gradients cause a significant jump in the image shift 
that exceeds the beam size. 

The above simulation demonstrates that the level of detail in 
the optically thin part of the jet is crucial for the reliability of the 
cross-correlation based image registration. If the jet has knots or 
bends, cross-correlation is rather robust against possible spectral 
index gradients. On the other hand, the method clearly does not 
work with smooth, straight jets that exhibit spectral index gra- 
dients. Also, the direction of the erroneous image shift due to 
spectral index gradient along these smooth, straight jets is such 
that it can mimic a true core shift. Therefore, we have repeated 
the core shift analysis in a "clean sample" from which such jets 
are removed. Statistics on the clean sample that comprises 94 
sources, however, did not show lower median values for 15.4- 
8.1, 15.4-8.4, and 15.4-12,1 GHz core shift distributions. This 
indicates that the effect of a possible overestimation of image 
shifts (and consequently core shifts) owing to a spectral index 
gradient along the jet is weak. 

We also estimated positional errors (from the imag e plane) of 
the co re components using the relation suggested by iFomalontl 
(119991) : (W ■ RMS)/(2P), where W and P are the convolved size 



and peak intensity of the component, RMS is the post-fit r.m.s. 
error associated with the pixels in a nine-beam area region un- 
der the component in the residual image. Because the cores are 
bright and compact, their signal-to-noise ratio S/N - P/RMS 
values are high, with a median value of a few hundred, making 
the corresponding positional errors small, with a median value 
of about 1 fias. This implies that the core shift error is typically 
dominated by the uncertainty of the image alignment for a sam- 
ple of core-dominated AGN jets. 

The typical accuracy level of about 50 //as achieved in our 
core shift measureme nts is comparable to th at from the self- 
referencing approach dSokolovskv et alj|201 lb . where the errors 
are dominated by position uncertainties of a referencing jet com- 
ponent, which is usually larger in size and weaker in flux den- 
sity with respect to the core. Phase-referencing observations may 
provide slightly better accuracy, down to 30 /uas, f or the 15.4 



8.4 GHz core shift, as reported bv lHada et al J ( 201lh . but cannot 
be used for a large number of sources as discussed in Sect. Q] 
Another complication of the phase-referencing method comes 
from the fact that the calibrator has its own core shift. 

5.3. Stationary jet feature problem 

Jets can exhibit standing features, like re-collimation shocks, in 
or close to their core region and these features could in princi- 
ple introduce a level of artificial core shift owing to the degra- 
dation of the angular resolution with wavelength. To study how 
strong an effect these standing features could have, we used the 
AIPS task UVMOD to create simulated VLBA data sets at 8.1, 

8.4, 15.3, 23.8, and 43.2 GHz. We used a real multi-wavelength 
VLBA observation to provide the (m, v) plane sampling and 
noise properties at every frequency, and as an input model for 
UVMOD we used three Gaussian components representing the 
core, a standing feature, and a jet feature. A set of simulated 
(u, v) data was generated with different flux ratios between the 
core and the standing feature (ranging from 10% to 50%) and 
with different distances between the two (0. 15 mas and 0.3 mas). 
As a result, we had a simulated multi-wavelength VLBA data 
set in which all wavelengths had exactly the same sky brightness 
distribution, but different (u, v) sampling and noise. Any appar- 
ent core-shift in this data set would be purely caused by the (u, v) 
coverage differences at different frequencies. Analyzing the sim- 
ulated data sets, we found that a close (within 0.3 mas from the 
43 GHz core) standing jet feature can contribute to the expected 
core shift effect between 43 GHz and 8 GHz owing to frequency- 
dependent blending on the level of ~10% of the expected core 
shift if the flux density ratio of the stationary feature and the 
core S S JS C ~ 10%, and reaching up to ~30% of the expected 
core shift if S st /S c ~ 50%. Between the frequencies used in 
this study, i.e., 15 and 8 GHz, the core shift is affected by the 
frequency-dependent blending by < 20 //as, which is less than 
the estimated accuracy of our measurements. We note that these 
simulations describe a very simple situation and a more detailed 
study of the effect of stationary features to the measured core 
shift using real data is warranted. However, this is beyond the 
scope of the current work. 

6. Jet physics with core shifts 

To increase the robustness of our subsequent calculations, from 
now on we use a vector average for each pair of the 15.4-8.1 and 
15.4-8.4 GHz core shifts. Moreover, if the core shift magnitude 
is smaller than 1 sigma level, we set Ar core i5_g = <x 15 _ 8 =51 //as 
to be used as the upper limit, since eliminating the small core 
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Figure 7. Distribution of redshift (left) and core shift measure 
Q rv (right) for 136 sources with the core shift derived between 15 
and 8 GHz. Gray filled bins represent 98 quasars (top), dashed 
bins represent 28 BL Lacs (middle), and cross-hatched bins rep- 
resent 10 galaxies (bottom). Empty bins show upper limits. 



shifts would introduce a bias. We excluded only the core shift 
vectors with |# cs - #j et | > 90°, which are likely to be dominated 
by errors. In total, we have core shifts between 15 and 8 GHz 
for 136 sources, 108 of which have both redshift (Fig. 13 left) 
and apparent jet speed measurements. If a source had a second 
epoch, we selected the one at which the dynamic range of the 1 5 
GHz image was higher. 

As shown by lLobanovl {1998!), core shift measurements can 
be used for deriving a variety of physical conditions in the 
compact jets. In particular, assuming equipartition between the 
particle and magnetic field energy density (k r = 1) and a jet 
spectral index aj et = -0.5 (S oc v +a ), the magnetic field in 
Gauss at 1 pc of actual distance from the jet ve rtex can be cal- 
culated through the followin g proportionality dHirotanil 120051 : 
|Q' Sullivan & Gabuzdall2009h 



Bi = 0.025 



n; v (i+z) 2 

<p5 2 sin 2 6» 



2\l/4 



(2) 




N 




BL Lacs 



B ipc (G) 



Figure 8. Distribution of magnetic field at a distance of 1 pc 
from the black hole for 84 quasars (top), 18 BL Lacs (bottom). 
Empty bins represent upper limits. 



where <p is the half jet opening angle, 9 is the viewing angle, 5 is 
the Doppler fac tor, and Q. rv is the core shift measure defined in 
lLobanovl d!998l) as 



Q„ = 4.85 • 10" 



where Ar c , 



Ar,, 



V\V 2 



(1+Z) 2 V 2 -Vl 



pc ■ GHz , 



(3) 



Vl is the core shift in milliarcseconds, and Di is 
the luminosity distance in parsecs. The calculated values of Q. rv 
in pc ■ GHz form a distribution ranging from 0.8 to 54.1 and 
peaking near the median of 13.6 (Fig. [7] right). The distribu- 
tions of £2 rv for quasars and BL Lacs are significantly differ- 
ent (p < 10 -4 ) as indicated by Gehan's general ized Wilcoxon 
test fr om the ASURV survival analysis package (lLavallev et al.l 
119921) . with medians of 18.6 and 7. 1 pc • GHz, respectively. 

The combination of ip, 5, and 9 in Eq. (0 is typically known 
for only a small fraction of sources, limiting the applicability 
of the formula. Therefore, the number of s ources in our sub- 
sample with known apparent jet speed B &w dLister et al.ll2009l) 
is larger by a factor of >2 than t he number of sourc es with 
known variability Doppler factor dHovatta et al.l 2009|). intrin- 
sic opening angle and viewing angle (e.g. lPushkarev et al.ll2.Q09l: 
ISavolainen et al.l2010l) . The denominator in (f2]i can be expressed 
through yS app only by substituting 5 = F~'(l - Bcos9)~ l and 
F m j n = ( 1 + 0l w ) x l 2 , and also taking into account that 2<p =; 
0.26 r 1 dPushkarev et al.l2009l) and B, ?? = Bsm9(\-Bcos9)- 1 . 
With these substitutions we are assuming the the jet is viewed at 
the critical angle 9 ^ F l that maximizes B app . We therefore ob- 
tain 



B\ * 0.042 Q^ 4 (l +z) 



1/2 



(4) 



where for /? app we use the fastest non-accelerating, r adial appar- 
ent speed measured in the source dListe r et al. 20091). 

First, we tested the consistency of B\ values calculated from 
Eqs. (f2]i and © for sources with previously measured 6 and 
Aipp values. For 40 sources out of 43 in common, comprising 35 
quasars and 8 BL Lacs, the results agree within the errors, with a 
median value of their ratio of 0.99. For three sources (0420-014, 
0804+499, and 1413+135), Eq. © gives several times higher 
values, most probably due to underestimated apparent speeds, 
to which (f2]i is more sensitive than d4), because the viewing an- 
gle 9 = arctan[2yS ap p(y6 2 pp + 6 2 - and the opening angle 
*f> — ^obs sin 9. Indeed, these sources have low apparent speeds 
but high Doppler factors, leading to low viewing angle estimates 
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Figure 9. Magnetic field at a distance of 1 pc from the central 
black hole versus fastest non-accelerating, radial apparent speed. 
The measurements are enveloped under the dashed aspect line. 
The middle dotted line shows the dependence based on the me- 
dians for redshift and core shift measure Q. rv . 

of 1°9, QP.2, 1°4 (ISavolainen et al.ll2010h . respectively, and in turn 
to overestimated magnetic field strengths. The uncertainties in 
B\ from Eq. [4] were calculated taking into account the errors in 
the core shifts, apparent speeds, and also from the assumption 
6 - F l , which is kn own to i ntrodu ce an additional source of 
errors, as discussed bv iListed (119991) . Since the aforementioned 
assumption is less correct for sources with low Lorentz factors, 
we excluded galaxies from the subsequent analysis. 

The distributions of B\ values derived from Eq. © for 84 
quasars and 18 BL Lacs shown in Fig.|8]have medians of 0.9^'^ 
and 0.4*2 j G, respectively, where the errors derived from a boot- 
strapping method indicate 95% confidence intervals. Gehan's 
generalized Wilcoxon test from the ASURV survival analysis 
package indicates that the distributions of B\ for quasars and 
BL Lacs are different at 99.6% confidence level. The difference 
is driven mainly by statistically higher redshifts for the quasars, 
and, though to a lesser degree, also by higher y6 ap p (Fig. |9). The 
B\ values would be comparable for these classes of objects, if the 
ratio of the median core shifts for BL Lacs to that of quasars was 
about 2.7, but that is not the case. Systematically stronger mag- 
netic fields in quasars can be a result of more massive black holes 
hosted by them and/or higher accretion rate in these objects, 
bacause larger black holes can accrete more matter, effectively 
powering the jets and accelerating particles to higher speeds. 
This scenario is also supported by an indication for quasars to 
have on average narrower intrinsic opening angle s than those of 
BL Lacs, as repor t ed by iPushkarev et al.l d2009l) . At the same 
time, IWoo&Urrvl (120021) argued that BL Lacs host black holes 
of comparable mass to quasars, though the estimates of black 
hole masses for quasars and BL Lacertae objects have been gen- 
erated by different methods, and this may introduce a bias while 
comparing the mass assessments. 

The absolute distance in parsecs of t he apparent VL B I core 
measured from the jet vertex is given bv (lLobanovl 19981) 

Q,. y Q,,(l + ^ pp ) 1/2 
fcoreW = — — ~ , (5) 

vsm8 v 

where v is the observed frequency in GHz. Then the correspond- 
ing magnetic field strength is B core = B\r^ K . In Fig. [TO] we plot 



Figure 10. Magnetic field in the 15 GHz core versus its depro- 
jected separation from the true jet base. The dashed line is the 
linear least-squares fit for quasar measurements. 

the derived r cole and B cole for the 15 GHz cores. Note that the 
core magnetic fields show a reverse tendency compared to the 
B\ values, with a median of 0.10 G for BL Lacs and 0.07 G for 
quasars, because their apparent cores are at different separations 
from the true jet base, with medians of 4.0, and 13.2 pc, respec- 
tively. 

On the projected plane (i.e., a VLBI map), the median sep- 
aration of the 15 GHz core from the central black hole is about 
1 pc, corresponding to about 0.14 mas for a source at z - 1. 
We summarize the derived physical quantities, such as Q, v , B\, 
B C0K , and r core in Table 0] 

The mass of a central bla ck hole can be r elated to B\, such 
that M bh « 2.7 x 1O 9 BiM (lLobanovlll998T) . Its typical value 
is ~10 9 Mq for quasars and BL Lacs. B\ can also be used to 
estimate the magnetic field Z?bh near the central black hole, as- 
suming a B oc R l dependence, where R is the half-width of 
the jet. Thus, we have Z?bh = B[Ri/Rbh, where R\ and Rbh 
are the the half-widths of the jet at 1 pc from the jet vertex and 
near the black hole, respectively. To derive R\, we assume that 
the jet is conical at distances larger than 1 pc. Then, we have 
^i = ^core _ ('"core - 1) tan <p, where R mie is the half- width of the 
15 GHz core, r core is the 15 GHz core separation from the jet 
vertex (both measured in pc), and tp is the intrinsic half jet open- 
ing angle. Because the black hole masses for BL Lacs are poorly 
known due to their weak emission lines, we restrict our anal- 
ysis to q uasars only. We deriv ed the following median values: 
<p = 0% (IPushkarev et al.ll2009l) . r core = 13.2 pc, R core = 0.35 pc. 
Then Ri as 0.2 pc and Z?bh ~ 2 x 10 3 G, assuming the jet 
width near the ~10 9 Mq black hole to be on the order of its grav- 
itational radius. The derived assessment of the magnetic field 
near the black hole is consistent with the theoretical value calcu- 
lated from a model of a thin, magnetically driven accretion disk 
dField & Rogerslll993l) . 

7. Summary 

We have implemented a method for measuring the frequency- 
dependent shift in absolute position of the parsec-scale core and 
applied it to multi-frequency (8.1, 8.4, 12.1, and 15.4 GHz) 
VLBA observations of 191 sources performed during 2006 
within the MOJAVE program. The method is based on results 
from (i) image registration achieved by a two-dimensional cross- 
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Table 4. Jet parameters derived for the sources with significant core shifts between 8 and 15 GHz. 



Source 


Epoch 


Opt. class 


z 


Papp 


Af corC] 15-8 GHz 


a n , 


Si 


^core 


'"core, 15.4GHz 












(mas) 


(pc GHz) 


(G) 


(G) 


(pc) 


0003-066 


2006-07-07 


B 


0.347 


2.89 


<0.051 


<4.55 


<0.20 


>0.22 


< 0.91 


0106+013 


2006-07-07 


Q 


2.099 


26.50 


<0.051 


<7.86 


<0.79 


>0.06 


<13.56 


0119+115 


2006-06-15 


Q 


0.570 


17.10 


0.324 


37.75 


1.63 


0.04 


42.07 


0133+476 


2006-08-09 


Q 


0.859 


12.98 


0.099 


13.69 


0.77 


0.07 


11.60 


0149+218 


2006-02-12 


Q 


1.320 


18.55 


0.196 


29.64 


1.69 


0.05 


35.83 


0202+149 


2006-09-06 


Q 


0.405 


6.41 


0.113 


10.87 


0.48 


0.10 


4.59 


0212+735 


2006-07-07 


Q 


2.367 


7.64 


0.143 


21.24 


1.27 


0.12 


10.65 


0215+015 


2006-04-28 


Q 


1.715 


34.16 


0.111 


17.08 


1.41 


0.04 


37.98 


0224+671 


2006-10-06 


Q 


0.523 


11.63 


0.139 


15.47 


0.75 


0.06 


11.75 


0234+285 


2006-09-06 


Q 


1.207 


12.26 


0.239 


35.82 


1.71 


0.06 


28.67 



Notes. Apparent speed values /3 app are taken from lLister et al .1 d2009j) . Table|4]is published in its entirety in the electronic version of Astronomy &■ 
Astrophysics. A portion is shown here for guidance regarding its form and content. 



correlation technique and (ii) structure model fitting. It has 
proved to be very effective and provided the core shifts in 163 
sources (85%), with a median of 128 /ias between 15 and 8 GHz, 
and 88 yt/as between 15 and 12 GHz. Despite the moderate sep- 
aration of the observing frequencies, the derived core shifts are 
significant (>2<x) in about 55% of cases, given an estimated typ- 
ical uncertainty of 50 fias between 15 and 8 GHz, and 35 /ias 
between 15 and 12 GHz. The errors are dominated by uncertain- 
ties from the two-dimentional cross-correlation procedure, be- 
cause the relative positional uncertainties of the compact bright 
cores are at a level of a few microarc seconds. The significant 
core shift vectors are found to be preferentially aligned with the 
median jet direction, departing from it by less than 30° in more 
than 90% of cases. 

We used the measured core shifts for constraining magnetic 
field strengths and core sizes for 89 sources. The magnetic field 
at a distance of 1 pc from the jet injection point is found to be 
~0.9 G for quasars and ~0.4 G for BL Lacs. Extrapolating all 
the way back by assuming a B oc R~ l dependence, the magnetic 
field in the close vicinity of the black hole is about 2 x 10 3 G. 
The core sizes, i.e., the distances from the true jet base to its 
apparent origin at 15.4 GHz are statistically larger in quasars 
than in BL Lacs, with a median of 13.2 and 4.0 pc, respectively. 
At these distances, the magnetic field has a median of 0.07 G for 
quasars and 0. 1 G for BL Lacertae objects. 

Future multi-epoch and multi-frequency VLBI observa- 
tions (including phase-referencing) of a pre-selected sample of 
sources with prominent core shifts are needed to address the 
question of the core shift variability, which can be used not only 
for astrophysical studies but also for astrometric applications. 
Ideally, these observations should be performed during and after 
strong nuclear flares that can be detected in advance in higher 
energy domains, e.g., optical or gamma-ray, and cover a wide 
range of observing frequencies, extending down to 5 or 2 GHz, 
where synchrotron self-absorption is essential. 
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Table 1. Derived core shift vectors. 
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43.2 
23.2 
88.3 
75.9 
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-65.1 
86.9 
-176.2 
-3.3 
171.0 
159.0 
-97.0 
1.5 
91.5 
178.0 

-104.5 
-138.0 
-12.7 
156.0 
173.1 
114.9 
118.5 
-3.1 
-149.5 
-105.3 
-94.9 
-87.1 
64.8 
-3.5 
-104.7 
15.8 
-175.3 
-174.0 
-124.9 
-111.2 
130.5 
131.8 
118.9 
140.8 
79.0 
70.4 
-65.2 
-20.7 
177.3 
147.2 
-47.7 
-79.1 
-55.0 
177.9 
-66.9 
135.0 
-55.7 



0.005 
0.018 
0.147 
0.201 
0.089 
0.168 
0.027 
0.074 
0.012 
0.040 
0.061 
0.043 
0.074 
0.218 
0.074 
0.147 
0.052 
0.209 

0.117 
0.141 
0.088 
0.174 
0.054 
0.078 
0.062 
0.074 
0.118 
0.015 
0.098 
0.118 
0.064 
0.034 
0.007 
0.114 
0.167 
0.166 
0.101 
0.238 
0.142 
0.047 
0.031 
0.126 
0.019 
0.007 
0.047 
0.190 
0.135 
0.178 
0.116 
0.006 
0.032 
0.094 
0.157 
0.013 
0.155 



-89.0 0.121 



160.4 
-45.3 
125.4 



0.200 
0.046 
0.062 



20.8 0.112 



-0.000 
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0.165 
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0.011 
0.016 
0.061 
0.043 
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0.217 
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0.132 
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0.112 
0.139 
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0.078 
0.061 
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0.014 
0.080 
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0.053 
0.026 
0.003 
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0.155 
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0.238 
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0.022 
0.120 

-0.009 
0.007 
0.040 
0.188 
0.134 
0.177 
0.114 

-0.002 
0.032 
0.088 
0.156 
0.005 
0.154 

0.121 

0.196 
0.046 
0.062 

0.109 
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Columns are as follows: (1) IAU name (B 1950.0); (2) epoch of observations; (3) 15.4 GHz median jet position angle; (4) position angle of the 
15.4-8.1 GHz core shift vector; (5) magnitude of the 15.4-8.1 GHz core shift vector; (6) 15.4—8.1 GHz core shift vector in projection on the 
median position angle; (7), (8), (9), and (10), (11), (12) the same as (4), (5), (6) but for 15.4-8.1 GHz and 15.4-12.1 GHz core shifts, respectively. 
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Table 4. Derived jet parameters. 
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